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Abstract 

This paper gives a momentum-space representation of the Argonne V18 potential as an expansion 
in products of spin-isospin operators with scalar coefficient functions of the momentum transfer. 
Two representations of the scalar coefficient functions for the strong part of the interaction are 
given. One is as an expansion in an orthonormal basis of rational functions and the other as 
an expansion in Chebyshev polynomials on different intervals. Both provide practical and efficient 
representations for computing the momentum-space potential that do not require integration or in- 
terpolation. Programs based on both expansions are available as supplementary material. Analytic 
expressions are given for the scalar coefficient functions of the Fourier transform of the electro- 
magnetic part of the Argonne V18. A simple method for computing the partial-wave projections 
of these interactions from the operator expressions is also given. 

PACS numbers: 21.45.Bc ,21.30.Cb 



I. INTRODUCTION 



The Argonne V18 potential {l| is one of a number of nucleon-nucleon interactions Qjljj^ 
that provide a quantitative description of experimental two-body observables below the pion- 
production threshold. It is distinguished from the other realistic interactions because it is 
expressed as an operator expansion with local configuration-space coefficient functions. This 
representation has advantages when used in variational Monte Carlo calculations. On the 
other hand, there are a number of calculations that require a realistic interaction that are 
more naturally performed in momentum space. These include some Faddeev calculations, 
relativistic few-body calculations, and calculations involving electromagnetic probes. In the 
momentum representation the variable conjugate to the relative coordinate is the momentum 
transfer. In calculations, both momenta appear, which requires either an interpolation or 
a separate Fourier transform for each pair of momenta. Fourier transforms of the V18 



potential have been used in some applications {4]. The purpose of this paper is to provide 
useful, tested and reproducible analytic approximations of the Fourier transform of the 
Argonne VI 8 potential for use in momentum-space calculations. The analytic forms allow 
for a direct calculation of the momentum-space interaction for any pair of initial and final 
momenta. In keeping with the traditional Argonne form, the momentum-space potential 
is given as a linear combination of products of spin-isospin operators with scalar functions 
of the momentum transfer. The resulting momentum-space potential has 24 terms. The 
additional six operators appear because the Fourier transform of the terms involving the 
operators L^V^(r) and (L ■ S)^Vi(r) each become a sum of two different momentum-space 
operators with different coefficient functions. In this work the Fourier transform is given 
for the strong part of the Argonne VI 8 potential, without the electromagnetic terms. This 
part of the potential must be treated numerically. The electromagnetic terms have analytic 
Fourier transforms, which are discussed in Appendix 3. The partial-wave projection of the 
momentum space potential is discussed in Appendix 2. It is constructed from the operator 
expressions by integrating over the angle between the initial and final momentum vectors, 
however unlike the configuration-space partial-wave projection, the integrals involve both 
the operator and the scalar coefficient functions. 



The Argonne VI 8 potential has the form 

18 

^ = 5ZK(r)0„ (1.1) 

n=l 

where Vn{r) are rotationally-invariant coefficient functions of the relative coordinate of the 
nucleons and the On are the eighteen spin-isospin operators given in Table 1., 



Table 1: Argonne VI 8 spin-isospin operators 
in coordinate-space 



Term 


snin-isosnin Operator in 

kjj XXX X ij k_J XXX V, y L/ V_jX CAj U V^X XX J 


r-snace 


Oi 


T 




yy2 








( fT^ ■ /Tn") 

\y 1 ^2)y 






\f I U 2)\i I 12) 




Or 


,S'n = .Sf/Ti ■ rlf^Tn ■ — 

'-'12 ^V'^l / 2 ''J 


^T-i ■ /To 


Oa 


^12 y 1 ' 2)y 




O7 






Oo 


yij iDjyTi T2) 






(L-L) 




Oio 


(L-L)(ri-T2) 




On 


(L ■L)(ai -(Ta) 




Ou 


(L ■ L)((Ti ■(T2)(ti ■ 


T2) 


On 


(L ■ Sf 




Ou 


(L-S)2(ri-T2) 






T12 = (3ri^r2^ - T ■ 






(ai ■ <T2)Ti2 






*S'l2ri2 




O18 


(riz + r2^) 





In this table T12 is the isotensor operator T12 := 3ri2r22 — Ti ■ T2. While the isospin 
operators, Tj, factor out of the Fourier transforms, the operators L^, L ■ S, (L ■ S)^ and the 
tensor operator 512 contribute to the Fourier transform. 

The Fourier transform of this potential can be expressed as a linear combination of 
24 momentum-space operators with scalar coefficient functions of the momentum transfer. 

3 



There are 24 operators because the L-L and (L-S)^ operators have two distinct contributions 
in momentum space. In appendix 1 it is shown that the potential matrix element (k'|\^|k), 
with q := k' — k, has the following five types of contributions: 



1. 



(2^ 



1 /" 1 /"°° 

_ j e-^iy-^yrv^^r)IdT = I— jo{qr)Vjirydr. (1.2) 



2. L-S 



(27r)3 
3. L L 



1 P 1 POO 

e-'^' ''Vj{r)L ■ Se'^'^dr = i(k x k') ■ S— - / ji{qr)Vj{rydr. (1.3) 

27r^g Jo 



^ J e-'^'-V,(r)L-Le^^-dr = 

1 /'CxD 1 poo 

-(k'xk)-(k'xk)-^ / j2 (gr )y,(r)r"dr + 2 (k'-k)-^ / j,iqr)V,irydr. (1.4) 
271" Q' Jo 2,71 q Jq 



4. (L ■ S)2 

/ e-^'^'-l-,(r)(L-S)V'^-dr = 

1 POO 1 /"CXD 

-(S-(kxk'))'-^ / j2(gr)\/,(r)r^dr + (k'xS)-(kxS)-^ / Mqr)V,{rydr. 

(1.5) 

5. S'i2 = 3(r ■ cri)(r ■ 0-2) - 0"i ■ 0^2 

^ 1 e-^'^'-'^Vir) (3(r ■ ai)(r • a,) - a, ■ a,) e^^ '^dr = 
-(3(q-ai)(q-a2)-gVi-a2)-^ / j2{qr)V{rydr. (1.6) 

These expressions are used to represent the momentum-space interaction as a sum of scalar 
functions of q := |q| multiplied by spin-isospin operators. These scalar coefficient functions 
of the momentum transfer that multiply the spin-isospin operators have the form of one of 
the integrals listed in Table 2: 



Table 2: Momentum-space scalar coefficient functions 



Scalar coefficient function 


dim 


indices 


VrrM ■= 2^/o°°Jo(gr)Kn(r)r2rfr 
■= 2^Io ji{qr)V„^{rydr 
Vm{q) := 2^ Jo j2{qr)Vm{r)r^dr 
Vm{q) := 2^52 /o j2{qr)Vm{rydr 


MeV fm3 
MeV fm^ 
MeV fm^ 
MeV fm^ 


m e {1,2,3,4,15,16,18} 
m e {7, 8, 96, 106, 116, 126, 136, 146} 
m G {9a, 10a, 11a, 12a, 13a, 14a} 
m e {5,6,17} 



where Vm{r) is the m*^ potential in the expansion (11. ip and Vma{q) and Vmb{q) are the two 
different functions that appear in fll.4p and (11.50 . These functions have finite limits as g — )■ 
in spite of the 1/g' coefficients since the Bessel function ji{qr) vanishes like as g — 0. The 
strong interaction contribution to the 24 scalar coefficients listed in Table 2 are numerically 
computed. The computational methods are discussed in section 3. Programs that compute 
these scalar coefficients are available as supplementary material to the electronic version 
of this paper. Quantities, like the binding energies in the test calculations, exhibit small 
sensitivities (in the sixth significant figure) to the precision of input constants. In the 
supplementary programs these constants are taken from the original VI 8 potential. 

The electromagnetic contribution to each of these operators can be represented in terms 
of known special functions. These contributions are important for precise low-energy cal- 
culations and can be added to the strong interaction coefficient functions when they are 
needed. The analytic expressions for the electromagnetic terms are given in Appendix 2. 

The resulting momentum-space potential has an operator expansion of the form 

(k'|i^|k) = J2^Uq)0m (1.7) 

where S = {1, 2, 3, 4, 5, 6, 7, 8, 9a, 96, 10a, 106, 11a, 116, 12a, 126, 13a, 136, 14a, 146, 15, 16, 17,18} 
and the 24 operators 0^ are given in Table 3. 
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Table 3: Argonne V18 momentum-space 



spin-isospin operators 



term 


spin-isospin operator 


Oi 


I 




O2 


(ti ■ T2) 




O3 


((Tl ■ 0-2) 




O4 


(CTi ■a2)(Ti ■ T2) 




O5 


- (3(q-<Ti)(q-0-2) - gVi ■ 




Oe 


- (3(q-(ri)(q-a2) - gVi ■ 


0-2) (ti ■ T2) 


O7 


i(k X k') ■ S 




Os 


^(k X k') ■ S(ri ■ T2) 






-(k' X k) ■ (k' X k) 




O96 


2(k' ■ k) 




OlQa 


-(k'xk)-(k'xk)(Ti-T2 


) 




2(k'-k)(Ti-T2) 






-(k' X k) ■ (k' X k)((Ti ■(T2 


) 


Onb 


2(k'-k)(ai-a2) 






-(k'xk)-(k'xk)(ai-a2)(Ti-T2) 




2(k'-k)((Ti-a2)(ri-T2) 






-(S- (kx k'))2 




Orsb 


(k' X S) ■ (k X S) 




Ol4a 


-(S ■ (k X k'))2(Ti ■ T2) 




Oub 


(k'xS)-(kxS)(Ti-T2) 




Oi5 


T12 




O16 


((Tl ■Cr2)Ti2 




Oi7 


-(3(q-ai)(q-o-2) -gVi ■ 


(^2) T12 


O18 


+ T2z). 





The Argonne V18 potential in momentum-space has the dimension MeV fm? . Dividing 
by he in Mev-fermi can be used to convert the momentum-space potential to a consistent 
set of units, (/m)^. 
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II. NUMERICAL FOURIER BESSEL TRANSFORMS 



This section summarizes an accurate numerical computation of the integrals in Table 2. 
These computations are used to test the accuracy of the approximations discussed in the 
next section. 

Because the configuration-space potential falls off asymptotically like e~"^^^, the radial 
integrals are evaluated with a finite cutoff at 20 fm. The Fourier-Bessel transforms are 
evaluated for momentum transfers q < 100 fm~^. With these cutoffs the maximum value of 
X := qr that can appear in the argument of the spherical Bessel functions in the integrals 
in Table 2. is Xmax = 2000. To evaluate these integrals the zeros of the spherical Bessel 
functions jo{x), ji{x), and j2{x) for < a; < 2000 are computed for each fixed value of q. 
For each value of q the integrals are expressed as sums of integrals between successive zeros 
of the spherical Bessel function that appear in the integral. If q is such that qr is never 
a zero of ji{qr) for < r < 20/m then the integral over r is performed using a 100 point 
Gauss- Legendre quadrature on the interval [0, 20/m]. If q is such that qr has zeros of ji{qr) 
for < r < 20 fm, then the integrals between zeros [gr^, grj+i] are computed using 20 Gauss- 
Legendre points when r^+i < 5 fm, 40 Gauss- Legendre points when 5 fm < rj+i < 10 fm 
and 80 Gauss- Legendre points when 10 fm < rj+i < 20 fm. For further details see~]5| 



III. APPROXIMATIONS 

This section discusses two approximations of the potential functions Vm{q) in Table 2 by 
expansions in known elementary functions. The first method approximates these potential 
functions by linear combinations of Chebyshev polynomials on three distinct intervals of 
momenta, for momenta up to 100/m~^. The second approach approximates these potential 
functions by a finite linear combination of orthonormal functions of the momentum transfer 
that have analytic Fourier-Bessel transforms. The configuration-space basis functions are 
associated Laguerre polynomials multiplied by decaying exponentials. These functions have 
analytic Fourier transforms that are rational functions of the momentum transfer In 
both approaches the coefficients of the expansion function are stored. The basis functions at 
any point can be generated efficiently by recursion and the potentials can be expressed as a 
finite linear combination of the basis functions. Both methods lead to efficient and accurate 



7 



approximations to the momentum-space potential. 

Figures 1 and 2 show the potential functions for the central and tensor parts (Vi(g) and 
V5(g)) of the interaction to illustrate the structure of typical potentials. 

A. Chebyshev expansions 

This section discusses the Chebyshev basis. The functions V„Jq) are replaced by a Cheby- 
shev polynomial approximation on the interval q E [a, 6] using |6| 

100 ^ ^ 2 

Vm{q)-Co/2 + y2cM-^ + - q) (3.1) 

^-^ — a — a 

n=l 

where 

Tn{x) = cos(n cos~"^(x))) (3.2) 

are Chebyshev polynomials and the coefficients Cn are computed using a Clenshaw-Curtiss 
quadrature js] : 

2 1 ~ -\-h h — 1 ~ 

Cn = j;^[^VUb) + ^-(^ + cos{rrj/N)) cos{njn/N) + (-)"-V;n(a)] (3.3) 

i=i 

with = 101. The functions Vm{q) are evaluated at the quadrature points qj := ^ -|- 
cos{nj /N) using the methods discussed above. This is repeated for q in each of three 
intervals, [a,b] = [0, 10], [10,50], [50, 100] and the 101 expansion coefficients associated with 
each of these three intervals are stored. The Chebyshev polynomials are computed using 
the recurrence relations 

Tn+i{x) = 2xTn{x) - Tn-i{x), Tq^x) = 1, Ti{x) = X. (3.4) 

For q larger than 100 fm~^ Vm{q) is approximated by 0. 

For the potentials V^iq), Vei^q) and Vn{q) it was necessary to add additional Cheby- 
shev expansions intervals between zero and ten /m~^. For V^lq) 21 polynomials were used 
on [0, .2]/m~^, 31 polynomials were used on [.2, .5]/m~^, 41 polynomials were used on 
[.5, 2.0]/m~^ and 71 polynomials were used on [2.0, 10.0]/m~^. For V^lq) 31 polynomials 
were used on [0, .5]/m~^, 41 polynomials were used on [.5, 2.0]/m~^ and 41 polynomials were 
used on [2.0, 10.0]/m~^. Similarly for Vi^iq) 31 polynomials were used on [0,1.0]/m~^, 51 
polynomials were used on [1.0, 5.0]/m"^ and 51 polynomials were used on [5.0, 10.0]/m~^. 
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This method provides an accurate and efficient representation for computing a momentum 
space V18 interaction. One of the supplementary programs (chebyshev-argonne.c) uses this 
method to compute the 24 coefficient functions in Table 2. 



B. Rational basis functions 

While the method of the previous section gives accurate results, a more straightforward 
approach is to represent the potential directly as an expansion in basis functions that have 
analytic Fourier transforms. In order to represent the potential, each of the scalar potentials 
Vm{(l), is approximated by an expansion in known basis functions. A method to compute 
both the expansion coefficients and a recursion formula to compute basis functions are given 
below. 

The functions Vm{r), rVm{r), and r'^Vm{r) that appear in the integrands of the integrals in 
Table 2 are expanded using an orthonormal set of radial functions that have analytic Fourier- 
Bessel transforms j?]. These functions are associated Laguerre polynomials multiplied by 
decaying exponentials in configuration space. Their Fourier-Bessel transforms have power- 
law fall of in momentum space. In addition, they vanish at the origin in a manner that can 
be used to explicitly cancel the factors 1/q and that appear in the definitions of Vm in 
Table 2. Both sets of basis functions can be generated efficiently using recursion relations. 
The cancellation of the factors 1/q and can be directly incorporated into the recursion 
that generates the momentum-space basis functions so the final expression for the potential 
does not require a special treatment for q near 0. 

The radial basis functions for different values of / are given below. The dimensionless 
parameter x := Ar is used in the basis functions, where A is a scale parameter that can be 
chosen to improve efficiency. The parameterization of the Argonne V18 interaction uses the 
value A = 7{fm)~^. The configuration-space basis functions are 




(3.5) 



where 




(3.6) 
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and the normalization coefficient is 



These functions satisfy the orthogonahty relations 



')niir)(f)miir)r^dr = 6mn- (S.i 







They have analytic Fourier-Bessel transforms given by 

'2 



(q) = \ - ji{qr)(j)ni{r)r dr. (3.9) 







For y = q/A the (pniil) can be expressed in terms of Jacobi polynomials: 

Mq) = l^, 2fiv+2 -^'"^''^^(^) (3-10) 

with normalization coefficient 

^ r(n + / + |)r(n + / + |) 

2(2n + 2/ + 3) n!r(n + 2/ + 3) ^ ' 

and 

rfa + n + l) / ^ \ r(a + /3 + n + m + l) 



p(«,/3)| 



a; 



- y - ir. (3.12) 

n ^„ 2™r a + m + l ^ ^ ^ ^ 



n!r(a + /3 + n + l)^^ I ; 2™r(a + m + 1) 
These functions satisfy the orthogonality relations 

(i>nliQ)(t>:nliq)q'^dq = 5mn- (3.13) 





These basis functions can be generated by using the recursion formulas for the associated 
Laguerre functions and Jacobi polynomials 

(n + l)L^+i(x) = (2n + a + 1 - x)L^^(x) - (n + a)L^,(x) (3.14) 

and 

2{n + l){n + a + f3 + l){2n + a + /3)Pi+f (x) = 
[(2n + a + /3 + l)(a2 _ ^ ^((2n + a + /3)(2n + a + /3 + l)(2n + a + /3 + 2)]P^"'^)(x) 

- 2(n + a)(n + /3)(2n + a + /3 + 2)Pi"f (x). (3.15) 
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These recursion relations can be modified to incorporate the normahzation constants (13. 7p 
and (13. lip directly into the recursion. The recursion for the normalized radial basis functions 
with (x = Ar) is given by: 



v/(2/+ 1)! 



Mr) = ^-^^i^Mr) (3.17) 



2n + l + 2[-2x (n-l)(n + l + 2/) ^ 

(pnlir) = —(pn-l,l{r) - \ / , O , OA <Pn-2M)- 3.18 

Vn + 1 + 2ly/n y n(n + 2 + 2/) 

Similarly, the normalized momentum-space basis functions with {y = q/A) are generated by 
the recursion: 

Mq) = 

111 1,1 



v/(2/ + 3)! VA3 ^1 . . . 2i±3 ^1 ...2;±i + 1)2^+2 



1 7/2 — 1 / 2/ -I- ^ 



(2r2 + 2/ + 3)n(n + 2/ + 2) 



(2n + 2/ + l)(n + / + |)(n + / + i) 
{2n + 21 + 1)(2/ + 2)(y2 + 1) + {2n + 21 + l){2n + 2l){2n + 21 + 2){y'^ - I) ~ 



2n{n + 21 + 2){2n + 2l){y^ - 1) 



{2n + 2/ + 3)(n - l)n{n + 1 + 2/)(n + 2/ + 2) 
{2n + 21- l){n + I + |)(n + / + \){n + I + \){n + I ~ \) 

{n){n + 21 + 2){2n + 21) ^ ^ ' 

Replacing 0oi(?) in P-l^P by 0oz(g) := ^Qi{si)l(i given by 

0o;(?) = 

' ' ' ' A-,-^i--, (3.22) 



^^/^2/ + 3)!^/A3 2^+3 ^/1 2f+l (l/2 + l)2«+ 
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to start the recursion in equations f l3.20l) - fl3.2ip generates (t>ni{(l) '■= 4>ni{(l) / which are well- 
behaved as g — 7- 0. Seventy expansion coefficients are used to construct the momentum-space 
potential for each value of m 

1 

Cnm = 7r^ (pnoir)Vm{r)r^dr m G {1, 2, 3, 4, 15, 16, 18} (3.23) 

Jo 

1 1"^ 

Cnm = 7r^ (t)ni{r)Vm{r)r''dr m G {7, 8, 96, 106, 116, 126, 136, 146} (3.24) 

Cnm = — 5- / <t'n2ij)ymij)r'^dr 171 G {9a, 10a, 11a, 12a, 13a, 14a} (3.25) 
2vr^ Jo 

1 f°° 

Cnm = 7r^ <Pn2{r)Vm{r)r^dr mG {5, 6, 17}. (3.26) 

The integrals are approximated using an 80 point Gauss- Legendre quadrature between 
and 10/m. The basis functions (f)ni{r) are generated using (13.161137181) . The scale parameter 
in the recursion for (j)ni{f) is taken as A = 7fm~^. 

The 70x24 expansion coefficients Cnm are stored. The momentum-space potential func- 
tions are then given by 

70 

Vm{q) = ^Cnm(f)nl{q) (3.27) 
n=l 

where the reduced expansion functions 4>ni{(l) '■= 4'ni{(l)/(i are generated recursively using 
(I320E22]). 

The full momentum-space potential in operator form is given by 

y=Y. ^™(^)^- (3-28) 

where Om are the 24 operators in Table 3 and q = y/k"^ + k'"^ — 2k' ■ k. 

One of the supplementary programs (rational-argonne.c) uses this method to compute 
the 24 coefficient functions in Table 2. 



IV. TESTS 

Two tests are performed on the potentials. First, the momentum-space coefficient func- 
tions, Vm{q), computed using the accurate numerical Fourier Bessel transforms, the Cheby- 
shev expansion and the rational basis function expansion are compared. For the second test 
both representations of potential are used to compute the deuteron binding energy and wave 
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functions. These results are compared to a direct calculation of these quantities using the 
partial-wave expansion of the original configuration space potential. 



The results of the first test are shown in Tables 4-7, which list values of the Fourier- 
Bessel transforms of the 24 radial functions computed using these three different methods 
for momentum transfers of 1,5,15 and 25 fm~^. 



These results are shown in Tables 4,5,6 and 7 for all 24 operators and a representative 
range of the momentum transfers. The columns labeled RFExp show the scalar potential 
functions using the rational function expansion, the columns labeled CExp show the same 
quantities using the Chebyshev expansion, while the columns labeled NET show the results 
of the direct numerical Eourier transform. Eigures 3-26 plot the difference of the approximate 
Eourier transforms with an accurate Eourier Bessel transform divided by half of the sum of 
these quantities. The solid curves are for the rational function expansion and the dotted 
curved are for the Chebyshev expansion. 
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Table 4: Values of scalar coefScients at 1 fm 
n RFExp CExp NFT 



1 


6, 


.789973e- 


-01 


6, 


.7899776- 


-01 


6, 


.7899776- 


-01 


2 


-4, 


.0193926- 


-01 


-4, 


.0193926- 


-01 


-4, 


.0193926- 


-01 


3 


-1, 


.692090e- 


-01 


-1, 


.6920906- 


-01 


-1, 


.6920906- 


-01 


4 


2, 


.358519e- 


-01 


2, 


.3567046- 


-01 


2, 


.3567056- 


-01 


5 


7, 


.2167396- 


-03 


7, 


.2182176- 


-03 


7, 


.2182176- 


-03 


6 


2, 


.8577326- 


-01 


2, 


.8604716- 


-01 


2, 


.8604676- 


-01 


7 


-5, 


.5115476- 


-01 


-5, 


.5115476- 


-01 


-5, 


.5115476- 


-01 


8 


-1, 


.6788886- 


-01 


-1, 


.6788886- 


-01 


-1, 


.6788886- 


-01 


9 


1, 


.7414156- 


-01 


1, 


.7414156- 


-01 


1, 


.7414156- 


-01 


10 


-3 


.2729886- 


-02 


-3 


.2729876- 


-02 


-3, 


.2729876- 


-02 


11 


1, 


.9991366- 


-02 


1, 


.9991366- 


-02 


1, 


.9991366- 


-02 


12 


-7, 


.4140606- 


-03 


-7, 


.4140606- 


-03 


-7, 


.4140606- 


-03 


13 


9, 


.0844226- 


-02 


9, 


.0844246- 


-02 


9, 


.0844246- 


-02 


14 


1, 


.2450176- 


-01 


1, 


.2450176- 


-01 


1, 


.2450176- 


-01 


15 


1, 


.1223886- 


-02 


1, 


.1223896- 


-02 


1, 


.1223896- 


-02 


16 


-1, 


.2149266- 


-02 


-1, 


.2160216- 


-02 


-1, 


.2160316- 


-02 


17 


2, 


.4032906- 


-03 


2, 


.4208186- 


-03 


2, 


.4208186- 


-03 


18 


6, 


.1249646- 


-03 


6, 


.1249646- 


-03 


6, 


.1249646- 


-03 


19 


1, 


.3042786- 


-02 


1, 


.3042746- 


-02 


1, 


.3042746- 


-02 


20 


-1, 


.7024096- 


-02 


-1, 


.7024016- 


-02 


-1, 


.7024016- 


-02 


21 


-7, 


.2272446- 


-03 


-7, 


.2272566- 


-03 


-7, 


.2272566- 


-03 


22 


-7, 


.8496866- 


-03 


-7, 


.8497076- 


-03 


-7, 


.8497076- 


-03 


23 


4, 


.5181936- 


-02 


4, 


.5182626- 


-02 


4, 


.5182626- 


-02 


24 


3, 


.9802516- 


-02 


3, 


.9802696- 


-02 


3, 


.9802696- 


-02 
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Table 5: Value of scalar coefficients at 5 fm 



n 




RFExp 




CExp 




NFT 


1 1.160699e+00 


1. 


,160699e- 


•00 


1, 


.160699e-00 


2 


-1, 


.360382e-02 


-1, 


.360382e- 


-02 


-1 


.360382e-02 


3 


-L 


.148807e-01 


-1, 


.148807e- 


-01 


-1 


.148807e-01 


4 


-1, 


.065288e-01 


-1, 


.065203e- 


-01 


-1, 


.065203e-01 


5 


4, 


.489757e-03 


A. 


.489763e- 


-03 


A. 


.489763e-03 


6 


4, 


.405849e-03 


4, 


.405371e- 


-03 


A. 


.405370e-03 


7 


-4, 


.623736e-02 


-A. 


.623736e- 


-02 


-A. 


.623736e-02 


8 


-1, 


.871380e-02 


-1, 


.871380e- 


-02 


-1 


.871380e-02 


9 


2, 


.471311e-02 


2, 


.471311e- 


-02 


2, 


.471311e-02 


10 


1, 


.480758e-03 


1, 


.480758e- 


-03 


1, 


.480758e-03 


11 


6, 


.027203e-03 


6, 


.0272036- 


-03 


6, 


.027203e-03 


12 


1, 


.465070e-03 


1, 


.465070e- 


-03 


1, 


.465070e-03 


13 


5, 


.222260e-03 


5, 


.222260e- 


-03 


5, 


.222260e-03 


14 


8, 


.233502e-03 


8, 


.233502e- 


-03 


8, 


.233502e-03 


15 


4, 


.828280e-03 


4, 


.828280e- 


-03 


A. 


.828280e-03 


16 


-4, 


.815794e-03 


-4. 


,815305e- 


•03 


-A. 


.815305e-03 


17 


1, 


.656921e-06 


1, 


.627533e- 


-06 


1. 


.627533e-06 


18 


4, 


.274306e-04 


4, 


.274306e- 


-04 


A. 


.274306e-04 


19 


4, 


.273833e-03 


A. 


.273832e- 


-03 


A. 


.273832e-03 


20 


1, 


.791462e-04 


1, 


.791461e- 


-04 


1. 


.791461e-04 


21 


9, 


.672551e-04 


9, 


.672550e- 


-04 


9, 


.672550e-04 


22 


1, 


.814761e-04 


1, 


.814761e- 


-04 


1, 


.814761e-04 


23 


1, 


.620319e-03 


1, 


.620318e- 


-03 


1, 


.620318e-03 


24 


1, 


.790086e-03 


1, 


.790085e- 


-03 


1, 


.790085e-03 
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Table 6: Value of scalar coefScients at 15 fm 
n RFExp CExp NFT 



1 


9, 


.321365e- 


-04 


9, 


.321031e- 


•04 


9, 


.321031e- 


•04 


2 


4, 


.123439e- 


-05 


4, 


.123387e- 


•05 


4, 


.123387e- 


•05 


3 


-1, 


.924812e- 


-05 


-1, 


.924669e- 


•05 


-1, 


.924669e- 


•05 


4 


-6, 


.648375e- 


-05 


-6, 


.643904e- 


•05 


-6, 


.643770e- 


•05 


5 


-9, 


.010902e- 


-06 


-9, 


.010512e- 


•06 


-9, 


.010512e- 


•06 


6 


1, 


.026393e- 


-05 


1, 


.026324e- 


•05 


1, 


.026323e- 


•05 


7 


5, 


.541260e- 


-06 


5, 


.540856e- 


•06 


5, 


.540856e- 


•06 


8 


2, 


.632043e- 


-06 


2, 


.631901e- 


•06 


2, 


.631901e- 


•06 


9 


-1, 


.962835e- 


-06 


-1, 


.962585e- 


•06 


-1, 


.962585e- 


•06 


10 


-9, 


.304609e- 


-07 


-9, 


.304846e- 


•07 


-9, 


.304846e- 


•07 


11 


-6, 


.015901e- 


-07 


-6, 


.015363e- 


•07 


-6, 


.015363e- 


•07 


12 


-1, 


.047669e- 


-07 


-1, 


.047529e- 


•07 


-1, 


.047529e- 


•07 


13 


-4, 


.725022e- 


-06 


-4, 


.725152e- 


•06 


-4, 


.725152e- 


•06 


14 


-1, 


.527634e- 


-06 


-1, 


.527584e- 


•06 


-1, 


.527584e- 


•06 


15 


2, 


.942747e- 


-06 


2, 


.942623e- 


•06 


2, 


.942623e- 


•06 


16 


-2, 


.895027e- 


-06 


-2, 


.892432e- 


•06 


-2, 


.892244e- 


•06 


17 


-2, 


.865458e- 


-10 


-3, 


.049671e- 


•10 


-3, 


.061292e- 


•10 


18 


9, 


.986465e- 


-08 


9, 


.985107e- 


•08 


9, 


.985107e- 


•08 


19 


-2, 


.604487e- 


-07 


-2, 


.604660e- 


•07 


-2, 


.604660e- 


•07 


20 


-6, 


.335039e- 


-08 


-6, 


.334951e- 


•08 


-6, 


.334951e- 


•08 


21 


-7, 


.055132e- 


-08 


-7, 


.055521e- 


•08 


-7, 


.055521e- 


•08 


22 


-1, 


.454468e- 


-08 


-1, 


.454569e- 


•08 


-1, 


.454569e- 


•08 


23 


-3 


.115148e- 


-07 


-3 


.115089e- 


•07 


-3 


.115089e- 


•07 


24 


-1, 


.394089e- 


-07 


-1, 


.394129e- 


•07 


-1, 


.394129e- 


•07 
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Table 7: Value of scalar coefScients at 25 fm 
n RFExp CExp NFT 



1 


-1, 


.386301e- 


•05 


-1, 


.383431e- 


•05 


-1, 


.383431e- 


•05 


2 


-6, 


.108349e- 


•08 


-6, 


.010007e- 


•08 


-6, 


.010007e- 


•08 


3 


8, 


.598072e- 


•07 


8, 


.595154e- 


•07 


8, 


.595154e- 


•07 


4 


1, 


.014189e- 


•06 


1, 


.003839e- 


•06 


1, 


.003915e- 


•06 


5 


-4, 


.600082e- 


•07 


-4, 


.599210e- 


•07 


-4, 


.599210e- 


•07 


6 


4, 


.739733e- 


•07 


4, 


.738710e- 


•07 


4, 


.73871 le- 


•07 


7 


2, 


.443040e- 


•08 


2, 


.442088e- 


•08 


2, 


.442088e- 


•08 


8 


9, 


.428095e- 


•09 


9, 


.412965e- 


•09 


9, 


.412965e- 


•09 


9 


-1, 


.534834e- 


•08 


-1, 


.533919e- 


•08 


-1, 


.533919e- 


•08 


10 


3, 


.457372e- 


•10 


3, 


.607579e- 


•10 


3, 


.607579e- 


•10 


11 


-3 


.619628e- 


•09 


-3 


.613201e- 


•09 


-3 


.613201e- 


•09 


12 


-1, 


.005784e- 


•09 


-1, 


.003137e- 


•09 


-L 


.003137e- 


•09 


13 


4, 


.666338e- 


•09 


4, 


.709390e- 


•09 


4, 


.709390e- 


•09 


14 


-3, 


.274714e- 


•09 


-3 


.270324e- 


•09 


-3, 


.270324e- 


•09 


15 


-5, 


.425469e- 


•08 


-5, 


.415686e- 


•08 


-5, 


.415686e- 


•08 


16 


5, 


.452722e- 


•08 


5, 


.398602e- 


•08 


5, 


.399643e- 


•08 


17 


-2, 


.888773e- 


•12 


-4, 


.263050e- 


•12 


-4, 


.213031e- 


•12 


18 


-5, 


.8521516- 


•09 


-5, 


.841440e- 


•09 


-5, 


.841440e- 


•09 


19 


-2, 


.512190e- 


•10 


-2, 


.555613e- 


•10 


-2, 


.555613e- 


•10 


20 


7, 


.827015e- 


•12 


7, 


.411785e- 


•12 


7, 


.411785e- 


•12 


21 


-5, 


.864134e- 


•11 


-5, 


.980761e- 


•11 


-5, 


.980761e- 


•11 


22 


-1, 


.617550e- 


•11 


-1, 


.652711e- 


•11 


-1, 


.652711e- 


•11 


23 


8, 


.297311e- 


•11 


8, 


.271942e- 


•11 


8, 


.271943e- 


•11 


24 


-5, 


.322500e- 


•11 


-5, 


.424546e- 


•11 


-5, 


.424546e- 


•11 



These tables show generally good agreement among the three methods of computation. 
At lfm~^ and 5fm~^ the Chebyshev expansion agrees with the direct numerical Fourier 
transform to between 5-7 significant figures for all 24 potentials. There is similar agreement 
at 15/m~^ and 25/m~^ except in potentials 17. The agreement between the potentials 
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calculated using the rational function expansion do not agree with the direct numerical 
Fourier transforms as well as the Chebyshev expansion. The accuracy depends on the 
particular potential and gets worse as the momentum transfer increases. Thus for precision 
calculations the Chebyshev expansion is preferred. 

Figures 1-24 provide a more complete picture of the nature of the errors in both approx- 
imations. Spikes in the errors occur near points where the potentials change sign. Some 
of the errors near zero are enhanced because the some of the plotted potential are divided 
by powers of the momentum transfer. For these terms the operators include compensating 
powers of the momentum transfer that vanish near the origin, so the contribution of the 
error in the full potential near the origin is reduced. The rational function expansions have 
larger relative errors near higher and lower values of the momentum transfer. This is not 
surprising because the basis functions are not local. The Chebyshev expansion is uniformly 
good, in part because it is a local expansion, so more intervals can be added as needed. The 
largest errors are in potential 17. At 10/m^^ its value is about —1.2 x 10~®, which is several 
orders of magnitude smaller than any of the other potentials at that momentum transfer. 



Table 8: deuteron s and d wave functions using Chebyshev expansion, 
rational function expansion and r-space partial waves 



kfm ^ 


Us{k)-CExp. 


Us{k)-RFExp 




Ud{k)-CExp, Ud(fc)-RFExp 


Ud{k)-pw 


0.0 


1.2695e+01 


1.2695e+01 1.2693e+01 


O.OOOOOe+00 


O.OOOOOe+00 0.00000e+ 


00 


0.5 


1.9609e+00 


1.9609e+00 


1.9609e+00 


-2.19827e-01 


-2.19811e-01 


-2.19826e- 


01 


1.0 


3.7684e-01 


3.7685e-01 


3.7684e-01 


-1.72164e-01 


-1.72159e-01 


-1.72164e- 


01 


1.5 


8.2472e-02 


8.2472e-02 


8.2471e-02 


-1.12429e-01 


-1.12429e-01 


-1.12429e- 


01 


2.0 


6.0809e-03 


6.0808e-03 


6.0806e-03 


-7.10857e-02 


-7.10863e-02 


-7.10857e- 


02 


2.5 


-1.3615e-02 


-1.3616e-02 


-1.3615e-02 


-4.45428e-02 


-4.45432e-02 


-4.45428e- 


02 


3.0 


-1.6153e-02 


-1.6153e-02 


-1.6153e-02 


-2.76853e-02 


-2.76854e-01 


-2.76853e- 


02 


3.5 


-1.3648e-02 


-1.3648e-02 


-1.3648e-02 


-1.69880e-02 


-1.69881e-02 


-1.69881e- 


02 


4.0 


-1.0153e-02 


-1.0153e-02 


-1.0153e-02 


-1.02233e-02 


-1.02234e-02 


-1.02234e- 


02 


4.5 


-6.9954e-03 


-6.9954e-03 


-6.9953e-03 


-5.98472e-03 


-5.98475e-03 


-5.98472e- 


03 


5.0 


-4.5270e-03 


-4.5270e-03 


-4.5270e-03 


-3.37040e-03 


-3.37043e-03 


-3.37041e- 


03 



As a second test the deuteron binding energy and wave functions are computed using 
the two different momentum-space potentials are compared to each other and to the di- 
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rect Fourier transform of wave functions computed using a configuration-space partial-wave 
calculation. 

The deuteron binding energy and the s and d wave functions are computed using the op- 
erator form of the Fourier transformed potential, by direct integration of the vector variables. 
The method of solution, which is discussed in |8|, uses the expansion (11. 7p directly with- 
out using partial waves. Calculations are performed for both the Chebyshev and rational- 
function representations of the momentum space potentials. 

These calculations are compared to a configuration-space partial-wave calculation. In 
that calculation, labeled pw in table 9, the wave functions are represented by an expansion 
in 70 configuration-space basis functions using the configuration-space basis functions (13.51) . 
Matrix elements of the partial wave projection of the Hamiltonian, with the configuration 
space Argonne VI 8 potential, are computed in this basis and the eigenvalue problem is 
solved numerically. The Fourier transform is computed by analytically Fourier transforming 
the basis functions. The solution of the eigenvalue problem gives an independent evaluation 
of both the binding energy and wave functions constructed directly from the configuration 
space potential. 

The deuteron binding energy obtained using the Chebyshev representation of the Fourier 
transform gives a deuteron binding energy of E^i = —2.242233 MeV. The rational function 
representation gives a d deuteron binding energy of = —2.242193 MeV compared with 
Ed = —2.242211 MeV using the configuration space partial-wave calculation. The binding 
energies based on all three calculation agree to within 22 eV. The computation used in 
the configuration-space partial-wave calculation is a Galerkin calculation and thus gives a 
variational bound on the binding energy. 

These eigenvalues differ from the experimental deuteron binding energy. This is because 
the momentum-space potentials used in these computations do not include electromagnetic 
corrections that appear in the Argonne VI 8 codes. The electromagnetic corrections con- 
tribute an additional +17.6 keV |9| to the binding energy which is consistent with the 
experimental binding energy of —2.2246 MeV. 

The s and d wave functions for all three calculations are compared in Table 10. The 
wave functions differ in the fifth or sixth significant figure, while binding energies of all three 
calculations differ in the sixth significant figure. 

The electromagnetic contributions to the Argonne V18 potential are important for low- 
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energy high-precision calculations. The Fourier transform of the electromagnetic contribu- 
tions of the Argonne V18 potential can be computed analytically, and can be added to the 
strong interaction contribution discussed in this paper when necessary. The analytic Fourier 
transform of the electromagnetic contribution is discussed in appendix 3. 

While this paper gives the momentum-space version of the operator expansion of the 
Argonne V18 potential, it is often useful to have partial-wave contributions of the potential. 
These can be computed from the operator matrix elements using a one-dimensional integra- 
tion over the cosine of the angle between the initial and final momenta. A simple method to 
compute the partial-wave projections from the operator expansion is given in the appendix 
2. 

The programs to compute the potentials Vm{<l) using both methods are available as 
supplementary material to the electronic version of this article. 

This work supported by the U.S. Department of Energy, contract # DE-FG02-86ER40286 
and National Science Foundation grants NSF-PHY-1005578 and NSF-PHY-1005501. The 
authors would like to acknowledge useful comments from Prof. C. Elster in preparing this 
manuscript. 
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V. APPENDIX 1 



In this appendix we compute the Fourier transform of the parts of the potential containing 
the five types of operators, I, L ■ S, L ■ L,(L ■ S)^, and S12 that appear in equations fll.3til.6p . 

L S: 

Let q := k' - k. 

The Fourier transform is 

^ J e-^^'-V,ir)S ■ (r x k)e^^-dr = j^, J e-^^-V,ir)S ■ (r x k)dr = 

A p 00 I 

7^ / E E i-^yM^^)yiU^)Yi*mmir)S ■ (r x k)dr. (5.1) 

) 1=0 m=-l 

Since r can be expanded as a hnear combination of spherical harmonics, Yim(r), the only 
terms that survive are the / = 1 terms. The integral over angles and the spherical harmonics 
replace f by q, giving 



POO 

/ Ji ■ (q X k.)r^dr 
Jo 



(2vr)3 Jo 

1 f°° 

zS ■ (k X k') X [—^ / ji{qr)Vj{rydr]. (5.2) 
27r q Jq 



Thus the Fourier transform has the form 



1 f e-*k'-ry^.(^)L ■ Se^^'^dr = zS ■ (k x k')/i(g) (5.3) 

J 



Hq) ■= ttV / Ji{qr)Vj{r)r'^dr. (5.4) 



(2vr) 

where 

hiq) ■■= . , 

The following relations are used to compute Fourier transforms of the remaining three 
operators: 

VJ(g) = f'{q)^ (5.5) 

Vj/(g) = fiq) + -/'(g) (5.6) 

q 

(a-V,)(b-V,)/(g) = 

^^^(r(g)-^) + ^/'(g). (5.7) 
q q q q 
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L L: 



The Fourier transform of this operator is 

^ I e-''^'-^V,{r){r x p) ■ (r x p)e^'^-^rfr = ^ / e-^'^'-^\/,(r)(r x k') ■ (r x k)e''^-^dr = 

(zV, X k') ■ (zV, X y V,{r)e-^^-dv = -(V, x k') ■ (V, x k)^ V,{r)j,{qr)rMr. 

(5.8) 

To compute the derivatives use 

(V, X k') • (V, X k) = (k . k')(V, ■ V,) - (k ■ V,)(k • V,) 

in the above to get 

- (V, X k') . (V, X k)^ / V,{r)Mqr)r'dr = - ((k' ■ k)Vj - (k' ■ V,)(k • V,)) /o(g) 

(5.9) 



(2vr) 

where 



An f°° 1 
Io{q) = J^:^ J Vj{r)jo{qr)r^dr = — / Vj{r)jo{qr)r^dr. (5.10) 



Evaluating this gives 



((k'-k)Vj-(k'.V,)(k.V,))/o(g) 



(k' ■ k)(/^(g) + ij^(g)) + (k^-qKk-q) (^.(^) _ 



To ehminate the derivatives use 



V{r){jQ{qr) - j'o{qr)—)r'^dr = ^ I V{r)j2{qr)r^dr = /2(g) 



and 



11/"°° 11 
m + -m = TT-, / V{r){f^{qr) - Uqr)- + 2Uqr)-)r'dr = 
q 2-71^ Jq qr qr 

V{r)j2{qr)r'^dr - — - / V{r)ji{qr)r^dr = ^{q) h{q). 



2vr2 7o ^^^^^^^ 2n^qJ, ' ''''^ ' q 

This gives 

-i^ / e-^'''-V,(r)(r x p) ■ (r x p)e*^-rfr = -(k' ■ k)(J2(g) - -h{q)) + ^^-^^^^^/2(g) 
(27r)^ J q q^ 

(5.12) 

which can be reexpressed in terms of cross products: 

^ ^ e-^'^'"y,(r)(rxp)-(rxp)e^'^"cir = -/2(g) ^ ' ^ + -(k^k)/i(g). (5.13) 



(27r)3y ^ ^^"^ g2 
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(L-S)^: 



The Fourier transform is 

^ / e"^('^'-'^)-r,(r)(L ■ S)'dr = J e-^^^'-^>^V,{r){S ■ (r x p)ydr 

Air f 



-((k'xS)-V,)((kxS)-V,)/o(g) = 
-((k' X S) • q((k X S) . q))^(/^(g) - ^/^(g)) - (k' X S) ■ (k X S)i/^(g) = 

- ((k' X S) ■ q((k X S) ■ q)) = ^hiq) + (k' X S) ■ (k X S)i/i(g) (5.14) 



which gives 



- ((S ■ (k X k')y\l2{q) + (k' X S) ■ (k X S)-Ji(g) (5.15) 



or 



^ I e-('''"'')y,(r)(L-S)2rfr = -((S-(kxkO)^^/2(g) + (k'xS)-(kxS)^/i(g). (5.16) 



(27r)^ 

1512 = (r ■(ri)(r ■0-2) - ■ 0-2): 

The Fourier transform is 

1 r / 1 



(^(V,-ai)(V,-a2)--(Ti-(r2V^J ^ y \/(r)-jo(gr)dr 



q q q q 

+ L,-a2{ll{q) + -lUq)) (5.17) 
o q 



Thus, 



1 /■ e-i''-l/(r) (3(r ■ a,){i- ■ a^) - a, ■ a^) e'^^ ^'dr = - -a.-a^] h-{q) 



(27r)3 7 v.vv X \ q q 

(5.18) 

where 

47r /" 

/o^g) = ^ j Vir)Uqrydr. (5.19) 
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VI. APPENDIX 2 



In this appendix we compute the partial-wave projection of potentials from the vector 
representation of the momentum-space potential. Using rotational invariance the partial- 
wave potentials can be expressed using a one-dimensional integral over the cosine of the 
angle between the initial and final momentum vectors. The method below is similar to the 



method first proposed in 



io|. 



Rotational invariance of the potential implies 

{j,fi,k,l,s\V\f,fi',k',l',s') 

= (j, /i, fc, /, s\U\R)VUiR) /i', k', I', s') 
= 5^ Di;{R)Di,^,{R){j, V, fc, /, s\V\2'. v\ k', s') 

vv' 

dRj2Di:,iR)Di,^,iR){j, u, k, I, s\V\j\ v\ k\ l\ s') 

vv' 

= 5^^.5,j,{k,l,s\\V^\\k',l\s') (6.1) 

where we have integrated over the SU(2) Haar measure with normalization J dR = 1 and 
defined the partial-wave potentials by 

{k,l,s\\Vmk',l',s'):=j^Yl E (6.2) 

IJ.=—j mim'j m.sm'g 

This kernel is rotationally invariant. Formally the partial- wave potential is 

{k,l,s\\Vmk'J',s') 

/I ^ \ \ 

x(k,/ii,/i2||V'||k',/i'i,/i2) 
X (^>/^'i7 ^>/^2l■5^"^l)(^^^>■5',m'Jj,/i)l^/^;(k'). (6.3) 
For any fixed rotation, R, rotational invariance of V gives 

(k,/ii,/x2||v^l|k',/.;,/.^) 
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\\U\R)VU{R)\\k',fi[,fi',) 

Using this expression in eq. (16. 3p gives 

{k,l,s\\Vmk'J',s') 

/I ^ \ \ 



(6.4) 



1 1 



,1/2 



(6.5) 

Next we eliminate three of the integrals in the potential matrix. For any fixed k' there is 
an R such that R{k')k' = z. Obviously R-^{k') = R,{(p')Ry{9')R^{^) where {9',(p') are the 
polar angles of k' and C, is arbitrary has this property where 



R. 



( cos{e) sin(e) ^ 

1 
\^ - sin(^) cos(^) J 

\ 



(6.6) 



In this case 



and 



For fixed k' define k" by 



cos(</)) — sin(0) 

sin(0) cos(0) 

V 

Rk' = k'z 



Rk = R;\ORy\0')R;\<p')k. 



k" = R-H0')R:H<l)')k. 



(6.7) 



(6.8) 



(6.9) 



(6.10) 



Given these fixed (primed) angles we change the unprimed integration variables k — k". 
We also have 



i?k = i?;i(e)k" 



(6.11) 
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We are also free to choose the angle ^ in R^^{^). We choose it so it transforms k" to the 
X — z plane. This is achieved by letting ^ be the azimuthal angle of k'" 

^ cos(0") sin(0") 0^ 
-UM>^ _ sin(0") cos(0") 
1 



R: 



( 



cos(0") sin(0") 



V 



(6.12) 



cos(0") 
y 



A;"sin(^") cos(0" 



- sin(0") 

\ ° 

With these substitutions the partial wave integral becomes 



k" sin(0") sin(0") 
\^ A;"cos(^"^ 



( k" sin(^") ^ 


\k"co^{Q") j 



(6.13) 



{k,l,s\\V^\\k\l\s') 
j dkW^-^ ^(^,^i,^,/X2k,m,)(/,m^.,m,|j- ^ 



xD*;/;^{R)D*;/;^{R){k'\ksin{e'') + zcos{^^^^^^ 

/I 1 1 

c^kW^— Y Yl < 2 ' 2 ' ""'^ '"^^-('^ sin(^")+z cos(^"))/^^»„^ (i?) 



X z^;;/; (i?)z};y;^ (i?) (a:"(x sin(^") + z cos(^") , fi'm\k% f^'i, /^2")^;('x 

X {l.f^[,l.f^2W,<){l',rn\,s',m',\jfi)Yi,^>^{^^ (6.14) 

Using properties of Clebsch-Gordan coefficients (i.e. D{R) < \ >=< \ > D{R) ® D{R)) we 
get 

{kJ,s\\Vmk'J',s') 



*l/2 



,1/2 



,1/2 



M=-J 



X (A;"(xsin(^") + z cos(^"), /ii, /iail V^H^'z, fx'^) 
X /i'l, ^, /i'sls', m',) (/', m;, s', m'Jj, (z)/^;...^(i?) 

/I ^ I X 

dk"dk' + 1 5Z ^2' 2' "^"'^ ^'^-^^ /^")^m,(^si'^(^") + ^ cos(6'")) 



M"=-j 
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X {k"{ksm{e") + z cos(^"), /ii, fi2\\V\\k'z, jj,'^) 

Since all of the dependence on </>', 6', (p" is in R, and the integrand is independent of R, these 
three angular integrals can be computed giving the multiplicative phase space factor of Stt^. 
What remains is the following integral over the cosine of the angle between the final and 
initial momentum: 

{k,l,s\\V^\\k',l',s') 

/I g^2 I 1 

du" . (9'/^!' 2'^2|s,m,)(/,m,,s,m,|j,/i")Y,,*„^(xv^ 

X (^>/^l.^./^2ls''"^s)(^''"^^s''Kb>'')^'m;(z)- (6-16) 

The last thing that needs for be addressed for an explicit formula is the spherical harmonics 



+ ZU") 



Y,U±Vl^-Hzu") = i-r^r-^f^p^Priu'') = f-^DURyiO")) (6.17) 



Yi>m[{z) = ^m'fi\j '^^ ■ (6.18) 
Using these in the above expression we are left with a 1 dimensional integral 

{k,l,s\\V^\\k\l\s') 



1 M"=-i 



X (A;"(xVl - u"^) + zm"), /ii, ^l2\\V\\k'z, /i^) 
1 1 I2I' + 1 



An 

Cleaning this up gives the following expression for the partial wave amplitude: 

{3,^^,k,l,s\V\3\^^\k\l\s') = 
5jy5,,>(k,ls\\V=\\k\l\s') 

with 

{k,l,s\\V^\\k\l\s') 
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(6.19) 



27 + lV 4:71 V 47r ^ " 2 

1 



X 



1 



du"Dl^o{Ry{cos-\u")){k"{±Vl - u"^) + zu"),fii,fi2\\V\\k'z,fi[,fi'^) 

X (^,/i'i,^,/i'2|s',m;)(/',0,s',ml|j/i") (6.20) 

where all repeated spin indicies are summed. This reduces the partial-wave integral to a 
one-dimensional integral. In this case there are no traces, and no momentum-dependent 
spin bases, but there are a number of spin sums. 

Explicit computation requires 

DLMcos-Hu")) ^ t ^7TlSElJfr^«ir-'fi;./^^"«&' (6.21) 

where all negative factorials are °° and 



R.,{cos-\u")) = I V ^ I (6.22) 



or 



VII. APPENDIX 3 

The electromagnetic corrections to the nucleon-nucleon interaction have the following 
forms for the pp, np, and nn systems: 

iVem,l{r) + VemA^) + VemA^) + VemA'^))^ + Vemfi{r)ai + t'em,9(r)5'i2 Vem,Vl{r)1^ " S (7.1) 
Vemnpir) = Vem,b{r)I + Vem,^{r)(Ji ■ (72 + Vem,ll{r)Si2 + tWl4(r)L ■ S (7.2) 
Vemnn{r) = Vem,7{r)cri ■ 02 + Vem,lo{r)Si2 + Vem,Vi{r)'h ■ S (7.3) 
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The Fourier transforms have the same structure as they do for the strong interactions. 

^emppiji) (^em,l 

Vemfi{q)^l ■ 02 + tW9(g)5'i2 + «fem,12(g)(k X k') ■ S. (7.4) 

Vemnp{q) = Vem^biq)! + Wem,8(g)cri " CF2 + Vem,ll{q)Sl2 + iiWem,14 (?) (k X k') ■ S. (7.5) 

(g)(kxk')-S. (7.6) 

where 

5*12 := -3((q ■ (Ti)(q ■ 0-2) - 0"i ■ 0-2)- (7.7) 
The coefficients of J,o"i ■ <r2 are 



00 



Wem,n(g) = 7T-^ / jo{qr)vem,n{r)r^dr n G {1, 2, 3, 4, 5, 6, 7, 8}; (7.8) 



Jo 

the coefficients of z(k x k') ■ S are 



1 f°° 

VemAq) = ^ J Ji{qr)vem,n{r)r'^dr n E {12, 13, 14} (7.9) 



and the coefficients of Su, are 



1 f°° 

VemAq) = 099 / j2{qr)Vem,n{^y^dr ne{9,10, 11} (7.10) 



where the individual terms are 



Vem,i{r) = ^(1 - (1 + ^br + ^{brf + ^(fer)3)e-^'-), (7.11) 

VemAr) = -^^(1 + + ^(&0^)e-^1, (7.12) 

VemA^) = ^fer^^(l)^ (7.13) 

nip 



VemAir) = ^(-7 ~ ^ + I Har)\)vem,i{r), (7.14) 



?,3 

^^en^,5(r) = ahcf3— {15 + 15br + 6{brf + {br)')e-'^), (7.15) 
C)o4 
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^ ^ ^^":(l + 6r + -(6rf)e-^'-), 



6m2 16 



(7.16) 



,(^) = _^(M!/i^(i + 5r + -(6r)2)e-^'-), 
^ 6m2 16^ 3^ ^ ^ ^ 



(7.17) 



Gmnrrip 16 



(7.18) 



adlc)^!^^ 1 1111 



4mp 



VemAr) = -^^^^^(1 - (1 + 6r + \{hrf + l(6r)3 + l(6r)^ + ^(6r)^)e-^0, (7.20) 



.e™,n(r) = -^^^^^(l-(l + &- + ^(&O^ + ^(&O^ + ^(&r)^ + Y^(H^)e-^0, (7-21) 



a(/j,c)^(4/ip - l)/ip/in 1 

r'^ 



7 



{l-{l + hr + -{hrY + -{hrY)e 



(7.22) 



(7.23) 



and 



VemMr) = -^M!/f!i - (1 + 6r + l(6r)^ + ^(&r)3)e-^0 



(7.24) 



where b = 4.27 and a = me/{hc), fip = 2.7928474, and /i„ = -1.9130427, /3 = .0189. 

The Fourier transform of most of the terms in the potential can be computed using direct 
integration, the identities 



(7.25) 



and the following relation Uj, with u = 1 + \ and /i = n + 1/2, gives the relation 
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g' (n + iy. ^,n + l + l n + l + 2 21 + 3 2/,2n 

7^ — , — n — ^^^^-r/b) 



(7.26) 



^n+l+l (21 + 1] 

which is vahd for n + I > — 1 . 

The only integral that can not be computed using these formulas involves the |ln(/i;r)| 
term that appears in VemA{f), which is an approximation to the vacuum polarization cor- 
rection to the pp interaction. The required integrals, which are evaluated in this appendix 
are: 



jo{qr)\ \n{ar)\rdr = — [7 + 3 ln(g/a) — 2ci(g/a)], 
^ 



(7.27) 



b'^ + q- 



jo{qr)\\n{ar)\re ^'dr 



:[btan-\l) -q-f + l Hq^a^ + b^a^)]- 



q2 52 



^2 2i 



(7.28) 



. Q POO 

jo(qr)\\n(ar)\r'^e~'^^dr = —— / jo(qr)\\n(ar)\re~'^dr, 

ob Jo 



(7.29) 



) q2 roo 

jo(qr)\\n(ar)\r^e~^^'dr = —— / jQ(qr)\\n(ar)\re~^'^dr, 

ob^ Jo 



(7.30) 



and 



^3 

jo{qr)\hi{ar)\r^e~^''dr = - — 
d¥ Jo 



jo{qr)\\n{ar)\re ^dr 



j2{qr)r ^dr = limji(x)/x = 1/3. 

x—^O 



(7.31) 



(7.32) 



The vacuum polarization integral only appears in the Coulomb potential which has the 
approximate form given in 12| - this approximation is adequate for binding energy calcula- 
tions. It appears in the following contribution to the proton-proton interaction 

^^en^(4)(r) = — (-7-g + |ln(ar)| + ^) — (1-e ^-{l + —br + — {brf + — {brf). (7.33) 



The Fourier Bessel transform of this interaction, 

1 



27r^g Jo 



ji){qr)Vem{i){r)r dr 
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1 



OTT O O iO iO 4o 

(7.34) 



can be computed analytically. There are three types of contributions 

"OO 



sin(gr)^(-7 - hc^hc{l - e-'^il + ^br + ^{brf + ^ibr)')dr, (7.35) 
OTT o lo lo 48 



(//) = / sin(gr) — |ln(ar)|a;ic(l-e-''^(l + -6r + -(6r)2 + -(6r)3),rfr, (7.36) 



and 
(///) 



27r2g2 



fivrA-r 11 3 1 

sin(gr) — — -a;ic(l - e^^Vl + — 6r + — (fer)^ + —(brf)dr. (7.37) 
^37r 8 ^ ^ 16 16^ ^ 48^ ^ ^ ^ ^ 



Integrals of the form (I) and (III) have the same form as the integrals discussed above. To 
calculate the integral (II) first replace r' = ar ^ r = r'/a to get 



1 /"^ 9n/ 

(//) = ^^J sm{qr'/a) — \ln{r')\ahc{l-e-'^'/''{l + 



^„,/„,^ 11 6r' 3 ,br' .0 1 /br' .... , , 
fr /"n J 1 — f — ) ))dr' 

16 a 16^ a ^ 48^ a ^ " 



1 2a 
2a7r2g2 37J- 



1 2a 



ahc I sm{qr' /a)\\n(r')\dr' 
d' 



11 d 3 r/2 1 ri'^ r°° 

(7.38) 



2a7r2g2 37J- 

Two integrals need to be performed to compute this term. They are 

sin(gr'/ a) | In(r') | dr' 



and 



(7.39) 



(7.40) 



roc 

/ sin(gr7a)e-^"'/"|ln(r')Mr'. 
Jo 

The integral 

POO 

/ sin(gr'/a)| ln(r')|(ir' = 
Jo 

— I sin(gr'/a) ln(r')(ir' + lim[ / sin(gr'/a)e~'^^ \n.{r')dr' — I sin(gr'/a) ln(r')(ir']. (7.41) 
Jo '^-^o ip-i Jo 

These integrals can be found in 13| : 

/ sin(gx) \n.{x)dx = — [7 + ln(g) — ci(g)] 
Jo Q 



where 



ci(x) = Ci(x) 



cos(t) 



(7.42) 



(7.43) 
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and 



f e-''"sin(gr)ln(r)c/r = — ^[6tan-^(|)-g7 + ^ln(g2 + 62)]_ (744) 
Jo 0^ + q b 2 

The quantity 7 is the Euler constant. 

Thus the first of the required integrals needed to compute the vacuum polarization con- 
tribution is 

POO 

I ln(r')| sm[qr' /a)dr' = 



/•l POO rl 

— / sin(gr'/a) ln(r')(ir' + lim[ / e"'^'' sin(gr'/a) ln(r')(ir' — / In(r') sin(gr'/a)] = 
Jo "'^o Jo Jo 

l*\ 1*00 

—2 / In(r') sin(gr'/a)cir' + lim / e""^^ sin(gr'/a) ln(r')cir' = 
Jo '^^0 Jo 

— [7 + ln(g/a) - ci(g/a)] + -[-7 + 1 hi{q^/a'')] = 
q q 2 

-[7 + 31n(g/a) - 2ci(g/a)]. (7.45) 

q 

Returning to the original expression - the first term in (II) is 

1 2a 



ahc / sm{qr' /a)\ln{r')\dr' 
Jo 

[7 + 3 ln(g/a) - 2ci(g/a)] (7.46) 



Sn'^q^ 

where 7 is the Euler constant. Or 

1^^^ sin(gr) ^ I In(ar) I a/ic = ^^j^ [7 + 3 ln(g/a) -2ci(g/a)]. (7.47) 
We also have to compute the second term in (II). The integral that is needed is 

sm(gr7a)| ln(r')|e ' [l + — h — ( — ) + 7^( — ))dr'. (7.48) 



37r^ag2 Jo 16 a 16 a 48 a 

This can be evaluated using 



00 

sin(gr7a)|ln(r')|e-^'"'/"dr' = 



00 roo 

sin(gr7a) ln(r')e-^'''/''dr' - 2 / sin(gr7a) ln(r')e-''"'/'^dr' (7.49) 



'0 Jl 

by differentiation with respect to h. The integral 



00 

sin(gr7a) \n{r')e-^''' /''dr' 
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2« Ji 



^(iq/a-b/ay _ ^{-iq / a-b / a)r' ^ ln(r')(ir' 



a /"°° dr' ^{-iq/a-b/a)r' 

+ 



2i Ji r' iq — b iq + b 



I b+iq ^ 
a > 



2i iq — b iq + b 

^ w T\rn/b-iq. .. .,^.b + iq. 
[{iq + b)Ei{ 1) + {iq - b)Ei{ ^] 



2i(g2 + 62) 



g2 + 62 I 



2i 



Thus 



sm{qr'/a) ln(r')e-^"'/'^rfr' 



g2 + 62 

Combining the two integrals gives 



2i 



sin(gr7a)|ln(r')|e-^'"'/"dr' 



g2 + 62 



(7.50) 



(7.51) 



(7.52) 



2 2i 

Putting all of the parts together gives the ln{r) contribution to the vacuum polarization 
contribution 



a^hc 
"irc^q^ 



[7 + 31n(g/a) -2ci(g/a)] + 



a^hc 



x{ 



g2 + fe2 



^[6tan-^(^) - g7 + I H^'/a' + b^la')] + 



2i 



(7.53) 
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This needs to be added to (I) and (II) to get the full vacuum polarization integral. These 
integrals can be computed using the methods used for all of the other potentials. 
The exponential integrals have simple derivatives 
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k[fm"'] 

FIG. 1. Vi{k). 
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Vg (k) 




FIG. 2. V5{k). 
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error (k) (rational-exact) vs (Chebyshev-exact) 
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FIG. 3. AVi rational (k), AVi Chebyshev{k) ■ 
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error (k) (rational-exact) vs (Chebyshev-exact) 




FIG. 4. AV2 rational {k), AV2 Chebyshev{k) ■ 
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error Vg (k) (rational-exact) vs (Chebyshev-exact) 
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FIG. 5. AV3 rational (k), AV3 Chebyshev{k) ■ 
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error V. (k) (rational-exact) vs (Chebyshev-exact) 



0.004 
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FIG. 6. AVirationalik), AViChebyshevik) ■ 
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error Vg (k) (rational-exact) vs (Chebyshev-exact) 




FIG. 7. AV5 rational (k), Chebyshev{k) ■ 
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error Vg (k) (rational-exact) vs (Chebyshev-exact) 




FIG. 8. AVio rational (k), AVq Chebyshev{k) ■ 
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error Vy (k) (rational-exact) vs (Chebyshev-exact) 
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FIG. 9. AVj rational (k), AV7 Chebyshev{k) ■ 
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error Vg (k) (rational-exact) vs (Chebyshev-exact) 
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FIG. 10. /^Vsrationalik), AVsChebyshev{k) . 
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error Vg (k) (rational-exact) vs (Chebyshev-exact) 




FIG. 11. AVghrationalik), ^VgbChebyshev{k) ■ 
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error V-|g (k) (rational-exact) vs (Chebyshev-exact) 




FIG. 12. AViobrational{k), ^ViobChebyshevik) ■ 
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error ^ (k) (rational-exact) vs (Chebyshev-exact) 




FIG. 13. AViibrational{k), ^VubChebyshevik) ■ 
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FIG. 14. AVi2brational{k), AVl2bChehyshev{k). 
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error V-|g (k) (rational-exact) vs (Chebyshev-exact) 
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FIG. 15. AVisbrational{k), ^Vi^bChebyshevik) ■ 
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error V^4 (k) (rational-exact) vs (Chebyshev-exact) 




FIG. 16. AVi4braUonal{k), ^VubChebyshevik) ■ 
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error V^g (k) (rational-exact) vs (Chebyshev-exact) 




FIG. 17. AVi5 rational{k), AV15 Cheby sheriff) ■ 
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error V^g (k) (rational-exact) vs (Chebyshev-exact) 




FIG. 18. AViQ rational{k), ^ViQChebyshevif^) ■ 
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error V^y (k) (rational-exact) vs (Chebyshev-exact) 




FIG. 19. AVij rational (k), AVn Chebyshevik) ■ 
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error V^g (k) (rational-exact) vs (Chebyshev-exact) 
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FIG. 20. AVi8 rational (k), ^VisChebyshevik) ■ 



55 



error V^g (k) (rational-exact) vs (Chebyshev-exact) 



- 




FIG. 21. AVgarational{k), AVgaChebyshev{k) ■ 
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error Vjq (k) (rational-exact) vs (Chebyshev-exact) 
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FIG. 22. AViOarational{k), AViQaChebyshev{k) ■ 
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error Vg^ (k) (rational-exact) vs (Chebyshev-exact) 




FIG. 23. AViiarational{k), ^VuaChebyshev{k) ■ 
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1e-05 - 



-1.5e-05 - 



error (k) (rational-exact) vs (Chebyshev-exact) 




FIG. 24. AVi2arational{k), AVl2aChebyshev{k) ■ 
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error Vgg (k) (rational-exact) vs (Chebyshev-exact) 



-0.0001 - 




FIG. 25. AViSarational{k), ^Vi3aChebyshev{k) ■ 
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error V24 (k) (rational-exact) vs (Chebyshev-exact) 




FIG. 26. AVi4arational{k), ^ViiaChebyshev{k) ■ 
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